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Abstract 

A measurement of the inclusive cross section for the process pp — > bbX — > ppX' at 
\fs = 7 TeV is presented, based on a data sample corresponding to an integrated lu- 
minosity of 27.9 pb~ collected by the CMS experiment at the LHC By selecting pairs 
of muons each with pseudorapidity \r\\ < 2.1, the value c(pp — >■ bbX — > jifiX') = 
26.4 ± 0.1 (stat.) ± 2.4 (syst.) ± 1.1 (lumi.) nb is obtained for muons with transverse 
momentum p T > 4GeV, and 5.12 ± 0.03 (stat.) ± 0.48 (syst.) ± 0.20 (lumi.)nb for 
px > 6GeV. These results are compared to QCD predictions at leading and next- 
to-leading orders. 
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1 Introduction 

The measurement of the cross section for inclusive b-quark production at the Large Hadron 
Collider (LHC) is a powerful probe of quantum chromodynamics (QCD) at very high energies. 
In addition, knowledge of the inclusive b-production rate from QCD processes helps to un- 
derstand the background in searches for massive particles decaying into b quarks, such as the 
Higgs boson or new heavy particles. 

The b-quark production cross section can be computed at next-to-leading order (NLO) in a 
perturbative QCD expansion HH31. The sizeable scale dependence of the result suggests that 
the contribution from the neglected higher-order terms is large [4—6] . The measurements per- 
formed at the Tevatron in pp collisions at y/s = 1.8 and 1.96 TeV [7. 8J, and at the LHC by 
the Compact Muon Solenoid (CMS) 19HTT) and LHCb fl2l [13| collaborations in pp collisions 
at \/s = 7 TeV in different rapidity ranges are generally consistent with the theoretical calcula- 
tions. However, the comparisons are affected by large theoretical uncertainties. 

The measurements of the cross section for the inclusive process pp — > bbX — >■ ppX' at y/s = 
7 TeV presented here allow for a comparison with QCD predictions in a kinematic domain 
where NLO calculations are more reliable because of the suppressed contribution of the gluon- 
splitting production mechanism (as discussed in [14] and the references therein). 

Experimentally, the dimuon final state allows for the selection of a sample with high bb event 
purity in the following wide kinematical region: muon pseudorapidity \tj\ < 2.1, where tj = 
— In [tan (9/2)] and 6 is the angle between the muon momentum and the counterclockwise 
beam direction, and muon momentum in the plane transverse to the beam axis px > 4 GeV or 
px > 6 GeV. Discrimination of the background from charm and light quark decays and from 
the Drell-Yan process is accomplished using the two-dimensional distribution of the two muon 
impact parameters {d xy ), defined as the distance of closest approach of each muon track to the 
interaction point projected onto the plane transverse to the beam axis. 

This paper is structured as follows. A brief description of the CMS detector is presented in 
Section |2] Section [3] describes the collision and simulated data used for this measurement and 
the selection criteria. Section [4] contains a detailed description of the categories in which events 
are grouped according to each muon's production process and kinematic features, while the 
fit to the impact parameter distributions is discussed in Section [5] Section [6] describes how 
the efficiency is computed and Section U\ is devoted to the determination of the systematic 
uncertainties. Section [H] reports the cross section measured in data and expected from QCD 
predictions. 

2 The CMS detector 

A detailed description of the CMS experiment can be found elsewhere [15J. The central fea- 
ture of the CMS apparatus is a superconducting 3.8 T solenoid of 6 m internal diameter. Within 
the field volume are the silicon tracker, the crystal electromagnetic calorimeter (ECAL), and 
the brass /scintillator hadron calorimeter (HCAL). Muons are detected in the pseudorapidity 
range \rj\ < 2 A by gaseous detectors utilizing three technologies: drift tubes (DT), cathode 
strip chambers (CSC), and resistive plate chambers (RPC), embedded in the steel return yoke. 
The silicon tracker is composed of pixel detectors (three barrel layers and two forward disks on 
either side of the detector, made of 66 million 100 pm x 150 pm pixels) followed by microstrip 
detectors (ten barrel layers, three inner disks and nine forward disks on either side of the de- 
tector, with the strip pitch between 80 and 180 pm). Thanks to the strong magnetic field and 
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high granularity of the silicon tracker, the transverse momentum pj of muons matched to re- 
constructed tracks is measured with the resolution better than 1.5% for pj < 100 GeV. The 
silicon tracker also provides the vertex position with ~15 f/m accuracy. The impact parameter 
resolution is measured with a sample of muons from Y(1S) — > }i + }i~ decays to be 28 pm and 
21 pm for muons with pj > 4 GeV and pj > 6 GeV, respectively. 

The first level (LI) of the CMS trigger system, composed of custom hardware processors, uses 
information from the calorimeters and muon detectors to select the most interesting events. 
The rapidity coverage of the LI muon triggers used in this analysis is\t]\ < 2.4. The high-level- 
trigger processor farm further decreases the event rate before data storage. 

3 Data selection and Monte Carlo simulation 

The data employed for this measurement were collected with the CMS detector during the 
2010 running period of the LHC. They correspond to an integrated luminosity £ = 27.9 ± 
1.1 pb~ [16|. A sample of events with two muons, each with transverse momentum pi > 
3 GeV were selected at the trigger level. Further requirements, designed to increase the purity 
of the muon candidates and to increase the fraction of muons from b decay in the sample, are 
applied at the analysis stage. A muon candidate is selected by matching information from the 
silicon tracker and muon chambers. The track must contain at least 12 hits from the silicon 
tracker, with signals in at least two pixel layers, and a normalized x 2 n °t exceeding 2. The 
overall x 2 obtained by combining the information from the tracker and the muon chambers 
should not exceed 10 times the number of degrees of freedom. Finally, each muon must be 
contained in the kinematical region defined by \rj\ < 2.1 and px > 4 GeV. We perform the 
measurement in this region and in a higher pj region where both the muons have pj > 6 GeV. 

Primary interaction vertices are reconstructed event-by-event from the reconstructed tracks. A 
candidate vertex is accepted if its fit has at least four degrees of freedom and its distance from 
the beam spot does not exceed 24 cm along the beam line and 1.8 cm in the plane transverse 
to the beams. Tracks are assigned to the primary vertex for which the track's distance to the 
vertex along the beam direction is smallest at the point of closest approach in the transverse 
plane. Muon tracks are required to have an impact parameter d xy perpendicular to the beam 
direction and with respect to its assigned primary vertex of less than 0.2 cm. Events are kept 
only if both muon tracks are assigned to the same primary vertex and both cross the beam axis 
within 1 cm of that vertex position along the beam direction. 

To remove muons from Z° decays, a selection on the dimuon mass M„„ < 70 GeV is applied. 
The mass range contributed by the Y resonances, 8.9 GeV < M /(f , < 10.6 GeV, is also rejected. 
Charmonium resonances and sequential semileptonic decays from a single b quark (for ex- 
ample b — > ]/ip X — > jifiX, or b — > cpX — > ppX') are rejected by removing dimuons with 
M w < 5 GeV. Events are selected if one and only one pair of muons is found satisfying all 
the criteria defined above. A total of 537734 events for pj > 4 GeV and 151314 events for 
pj > 6 GeV pass these requirements. 

Two samples of simulated Monte Carlo (MC) events were generated using the minimum-bias 
settings of PYTHIA 6.422 071 (parameter MSEL=1), with the Z2 tune flT8j H3, and incorpo- 
rating the CTEQ6L1 parton distribution functions (PDF) 11201 . To increase the generation ef- 
ficiency within the selected acceptance, a filter was applied at the generator level requiring 
two muons with p| en > 2.5 GeV and |// gen | < 2.5 for the measurement with pj > 4 GeV, or 
p T en > 5 GeV and |f/ gen | < 2.5 for the measurement with pj > 6 GeV. The generated samples 
include events with muons originating from the decay of light mesons (mostly charged pions 



and kaons) within the tracker volume. A third MC sample was produced to simulate the Drell- 
Yan process. MC events, including the full simulation of the CMS detector and trigger via the 
Geant4 package [21], are subjected to the same reconstruction and selection as the real data. 

4 Templates for different muon classes 

The fraction of signal events (pp — > bbX — > jijiX') in the data is obtained from a fit to the 2D 
distribution of the impact parameters of the two muons. For this purpose, reconstructed muons 
in the simulated events are separated into four different classes, defined according to their 
origin. The single-particle distributions of the transverse impact parameter d xy are obtained 
for each class from simulation and fit using analytical functions. From these functions, the 2D 
templates are built symmetrically. This procedure is described in the following section. 

4.1 Definition of muon classes 

Information from the generation process is used to assign each reconstructed muon in the simu- 
lation to a well-defined category. Reconstructed muon candidates are linked to the correspond- 
ing generated charged particles with a hit-based associator, which reduces the probability of 
incorrect associations to a negligible level. Tracks are assigned to one of the following classes: 

1. B-hadron decays (B): muons produced in the decay of a B hadron, including both direct 
decays (b -4 ]i~X) and cascade decays (b ->■ cX -> pX', b -4 tX ->■ jiX', b -4 }/\p X -4 
p±X); 

2. Charmed hadron decays (C): muons from the semileptonic decays of charmed hadrons 
produced promptly; 

3. Prompt tracks (P): candidates originating from the primary vertex, mostly muons from 
the Drell-Yan process and quarkonia decays. This category also includes punch through 
of primary hadrons, and muons from decays of charged pions and kaons in the volume 
between the silicon tracker and the muon chambers; 

4. Decays in flight (D): muons produced in decays of charged pions or kaons (which may 
come either from light- or heavy-flavor hadrons) in the silicon tracker volume. 

Table IT] gives the single-muon sample composition from the simulation for MC events passing 
the full selection and dimuon trigger. While the fraction of muons from decays in flight (D) 
decreases at larger pj, the prompt component (P) increases due to the Drell-Yan muons. 

The predicted composition of the dimuon events from the simulation is shown in Table|2l where 
PX is defined as the sum of the PB, PC, and PD contributions. The uncertainties given in the 
table are the statistical uncertainties from the simulated samples. 

Figure [l| shows the d X y distributions for muons with pj > 4GeV from the simulation for all the 
classes above except for the prompt tracks, where muons from decays of Y(1S) in the collision 
data are used after removing the background with a sideband subtraction technique. 

The prompt d xy distribution is fit with the sum of a Gaussian centered at zero and an exponen- 
tial function. This combination of functions accounts for the detector resolution effects. The 
distributions of the other classes are fit using, in addition, a second exponential term. The 
functions are shown by continuous black lines overlaid on the histograms in Fig. Ill while the 
black points represent the template histograms obtained by evaluating the fit functions at each 
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Table 1: Percentage of each muon class in the simulated events for two px requirements. The 



uncertainties are statistical only. 



Source 




Fraction in simulation (%) 






p T > 4GeV 


px > 6 GeV 


B hadron (B) 




77.8 ± 0.2 


79.8 ± 0.4 


C hadron (C) 




14.0 ±0.1 


12.6 ± 0.1 


Prompt sources 


(P) 


1.84 ± 0.04 


3.44 ± 0.08 


Decays in flight (D) 


6.37 ±0.07 


4.21 ± 0.09 



Table 2: Percentage of dimuon event sources in the simulation for two different px require- 
ments. PX represents the sum of the contributions from PB, PC, and PD. The uncertainties are 
statistical only. 



Source 


Fraction in simulation (%) 




p T > 4 GeV 


px > 6 GeV 


BB 


71.6 ±0.2 


74.6 ± 0.4 


CC 


9.24 ± 0.08 


8.67 ±0.14 


BC 


5.66 ± 0.07 


5.22 ±0.11 


PP 


1.84 ±0.04 


3.43 ± 0.08 


DD 


1.49 ±0.04 


0.73 ± 0.04 


BD 


6.01 ±0.07 


4.40 ± 0.10 


CD 


3.69 ± 0.05 


2.53 ± 0.08 


PX 


0.48 ± 0.02 


0.40 ± 0.03 



bin center. The ratio of the MC distribution to the fit values are shown in the lower plots of 
Fig. FT] The templates for muons with px > 6 GeV are obtained in a similar way. 

4.2 Two-dimensional template distributions 

In principle, the dimuon events could be split into sixteen different categories by combining the 
four classes defined above for each muon. In order to reduce the number of categories to ten, 
the d X y distributions are symmetrized (i.e., BC=CB, BD=DB, etc.) using a method originally de- 
veloped by the CDF collaboration [8J. The one-dimensional (ID) histograms, built as described 
above, are normalized to unity within the fit range < d xy < 0.2 cm. The symmetrized 2D 
template histogram for the events with a muon of class p and another of class c (p, (X = 1, . . . , 4 
according to the definition in Section [4T|| is then constructed as: 



ij 



^sj + stsf), 



(1) 



where S^ is the content of the i bin of the histogram describing the class p, and analogously 
for index j and class a. In this way, ten symmetric distributions are obtained. In practice, the 
few events from the PX category are neglected, thus reducing the number of significant classes 
to seven. 

The ID projections of the seven templates are shown in Fig.|2]for muons with px > 4 GeV. 



5 Measurement of the sample composition 

Consistent with the symmetric 2D templates, the data events are also randomized by taking 
the impact parameters of the two muons in each event, and filling the bin corresponding to 
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Figure 1: Comparison, for muons with px > 4GeV, between the template d xy histogram (red) 
and the fitted function (black) for muons coming from B hadrons (class B, top left), charmed 
hadrons (class C, top right), prompt tracks (class P, bottom left), and decays in flight (class D, 
bottom right). The templates for B, C, and D come from simulation. For the prompt tracks, 
the distribution is obtained from data. An enlargement of the prompt-track distribution for 
dxy > 0.05 cm is shown on a linear scale as an insert in the lower-left plot. For each template, 
the ratio of the d xy histogram to the fitted function is shown at the bottom. 



\dxy(]i\), d X y(]i2)] or to [d X y(ji2), d X y(}ii)] according to the outcome of a random number gener- 
ator. 

The fractions of the individual contributions to the observed distribution are determined with 
a binned maximum-likelihood fit. The fit minimizes the function: 



21n(L) = -2<^ £ [mjin(lij)-k 



w=i 




(2) 
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Figure 2: ID projections of the d xy templates used in the fit for muons with pj > A GeV, for the 
BB, CC, PP, DD categories (left) and the BC, BD, CD ones (right). 



where nn is the content of the data histogram in the bin (i, j), Ijj = YLk[fk- Tk,ij\> where T\ is the A: th 
template (k = 1, . . . ,7), and fa is the fit parameter expressing the fraction of events from the k th 
source. The fitted fractions are subject to the normalization condition Yjk=\ /fc = 1- To reduce the 
number of fit parameters and ease the fit convergence, the three parameters /bC//bd/ and /qd 
are constrained so that the ratios /bc/ '/be, /bd/ '/bb, and feu/ fee are compatible with the MC 
expectations within their statistical uncertainties. In Eq. pi, k' is the index of the constrained 
templates (BC, BD, CD), r^ is the ratio of the constrained fit fraction with respect to the reference 
fit fraction (for instance in the BC case r^c = /bcZ/bb), r$ c is the ratio of the constrained 
fraction and reference fraction in the simulation, and <7^ c its statistical uncertainty from the 
number of simulated events. 

The BC component originates from the production of an extra cc pair from gluon splitting in a 
bb event. The production rate of cc pairs from gluon splitting has been measured at LEP [22- 
24], and found to be 50% higher than theoretical predictions |25|. The measured bb rate H26ti28| 
is about 10 times smaller and has a negligible effect on the BC component. In contrast, the BD 
and CD contributions are related to the misidentified muon rate in events with true B and C 
production. These rates are determined from the MC simulation, and have been checked using 
direct measurements in the data [29] . The systematic uncertainties on the fit constraints are 
discussed in Section 1731 

Table [3] gives the results of the fit to the data sample. The quoted uncertainties are obtained 
from the fit and are statistical only. The measured BB fraction is smaller than expected from the 
simulation, while the DD fraction is larger. Projections of the d xlJ distributions with the results 
of the fits are shown in Fig.[3]for the two pj selections. 



6 Efficiency determination 

The total efficiency e is defined as the fraction of signal events produced within the acceptance 
(pi > 4GeV or pj > 6GeV, \r]\ < 2.1 for each muon) that are retained in the analysis. In the 
simulation, the values of eMC = (44.3 ± 0.1)% and (69.9 ± 0.1)% are computed for signal events 
with a pi threshold of 4 and 6 GeV, respectively. 



Table 3: Results of the likelihood fit to data for the percentage of each dimuon source with two 
different muon pj requirements. The BC, BD, and CD fractions are constrained to their ratios 
to BB and CC fractions as expected from the simulation. 



Source 


p T > 4GeV 


p T > 6 GeV 


BB 


66.8 ± 0.3 


70.2 ± 0.3 


CC 


9.2 ±0.6 


5.5 ±1.2 


BC 


5.2 ±0.1 


4.9 ±0.1 


PP 


1.7 ±0.3 


4.0 ± 0.4 


DD 


7.8 ± 1.1 


9.5 ±2.1 


BD 


5.6 ±0.1 


4.2 ±0.1 


CD 


3.7 ±0.9 


1.6 ±0.5 
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Figure 3: Top: The projected d xy distributions from data with the results of the fit for muons 
with px > 4 GeV (left) and pj > 6 GeV (right). The distribution from each dimuon source is 
shown by the histograms. Bottom: The pull distribution from the fit. 



To compare these values to efficiencies measured in data, the selection procedure is divided 
into three steps, each defined relative to events passing the previous one: 

1. muon selection ("MuSel"): events having at least two selected muons, each associated 
with a reconstructed vertex; 

2. event selection ("EvSel"): events passing the dimuon invariant mass requirements, with 
both muons belonging to the same vertex; 

3. trigger selection ("Trg"): events passing the trigger requirements. 



The efficiencies obtained by counting the signal events passing each step in the simulation are 
given in Table El 

The total efficiency can alternatively be expressed on an event-by-event basis by defining the 
efficiency e, to select the z' th signal event as e ; - = e^MuSei ' e z',EvSei ' £ i,Trg- The (pj,f]) distri- 
bution of the signal events and the efficiency e^EvSei can only be extracted from simulation. 
The efficiencies e/,MuSei an d £ i,Trg can De found as the products of the single-muon efficiencies, 
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S{ = s fll (pj,r]) ■ s^{jpi,y]), under the assumption that the single-muon efficiencies e }li only 
depend on the px and rj of the muon. This assumption is found to be compatible with the 
efficiencies determined in the simulated sample, within their statistical uncertainties. 

Table 4: Efficiencies (in percent) at each step of the analysis found from the simulation and from 
the data. The last column reports the overall efficiency obtained from the product of the three 
efficiencies shown. The event selection efficiency ^EvSei cannot be found with the data, so the 
MC simulation value is used. The bias and feed-through corrections described in the text are 
also included in the overall efficiency. Only statistical uncertainties are reported. 



Sample 


e MuSel 


^EvSel 


6 Trg 


e 


p T >4GeV,MC 
p T > 4GeV, data 


64.8 ±0.1 
69.5 ± 3.6 


78.0 ± 0.1 


87.7 ±0.1 
86.1 ±2.0 


44.3 ±0.1 
48.8 ± 2.9 


p T >6GeV,MC 
px > 6 GeV, data 


83.6 ±0.1 
87.0 ± 3.4 


90.1 ±0.1 


92.8 ±0.1 
93.4 ±2.1 


69.9 ±0.1 
74.4 ± 3.8 



In the data, the single-muon selection and trigger efficiencies are measured in intervals of px 
and rj with the "tag-and-probe" (T&P) method H29j|30||, which employs a sample of ]/ip —■ p + p~ 
events selected with minimal trigger requirements. The selection efficiency found from this 
method is consistent with the value from the simulation (^MuSel^MuSel = 1-073 ± 0.054 for px > 
4 GeV and 1.041 ± 0.047 for p T > 6 GeV), as is the trigger efficiency (e^/e^g = 0-982 ± 0.028 
for p T > 4 GeV and 1.006 ± 0.023 for p T > 6 GeV). 

Differences in the kinematic distributions between the ]/tp sample and the bb events might 
imply different bin-averaged efficiencies, causing biases in the region close to the acceptance 
thresholds. An overall bias correction of 0.966 ± 0.015 (1.004 ± 0.012) is computed when com- 
paring the efficiencies in the simulation computed with the T&P method and those obtained 
with the signal in the range px > 4 GeV (px > 6 GeV), where the uncertainties are statistical 
only. 

Another correction to the total efficiency is applied to take into account the feed-through of 
events where one of the muons has a true px below the selection limit, whereas the recon- 
structed px is above it. This effect is computed using the simulation by finding the fraction 
of selected events with at least one muon generated outside of the acceptance, and is equal to 
0.990 (0.980) for p T > 4 (6) GeV, with negligible uncertainties. 

The overall efficiency is computed as the product of the efficiencies for the muon selection and 
trigger, as obtained with the T&P method in data, times the event selection efficiency found 
in the simulation, divided by the bias and the feed-through corrections. Results are shown in 
Table |U 



7 Systematic uncertainties 

Several sources of systematic uncertainties have been considered for this measurement. They 
are divided into uncertainties due to the model dependencies for both the signal and the back- 
grounds, the effects related to the impact parameter resolution, the fit method, and the mea- 
surement of the efficiency. Each of these is described separately in the subsections below. 

7.1 Model-dependent uncertainties 

The impact parameter projected onto the plane transverse to the beam axis of a muon produced 
in a hadron semileptonic decay is related to the parent hadron's proper decay time t by: 
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d X y = ^7 ct sir\5 sm9, (3) 

where {5 is the ratio of the hadron velocity and the speed of light c, 7 = (1 — /3 2 ) 1,/2 , 5 is the angle 
between the muon and the hadron directions in the laboratory frame, and 9 is the polar angle 
between the hadron direction and the beam axis. Uncertainties in the parent lifetime affect the 
simulation of the proper distance distribution ct, and uncertainties in the parent hadron energy 
spectrum affect the Lorentz boost factor fiy and the angle 5. The three general categories of 
systematic uncertainties due to these model dependencies are: 

b- and c-hadron properties: four of the long-lived B hadrons produced at the LHC decay to muons 
at a non-negligible rate. While the B^ and B u lifetimes are known with a precision better than 
1%, the B s and Aj, lifetimes are measured with larger uncertainties. Simulated MC events with 
B s and Ab decays are reweighted so as to vary the corresponding lifetimes by their uncertainties 
[l3"T| , the templates are recomputed, and the fit is repeated. The fit result changes by ±2.1% 
(±1.5%) for pi > 4 (6) GeV. The effects from uncertainties on the B^, B u , and c-hadron lifetimes, 
similarly evaluated, are negligible. The b-hadron sample composition has been measured by 
experiments at LEP and the Tevatron [31 J, and by LHCb |32]. While substantial agreement 
has been found for the B s / (B u + B^) fraction, a sizable discrepancy was observed for f^ b = 
Ab/ (B u + Bd). The results presented here are obtained using the averages between the LEP 
and LHCb results, f Ab = 0.18 ± 0.09 (0.165 ± 0.075) for p T > 4 (6) GeV, where the uncertainties 
correspond to half the difference between LEP and LHCb. Varying f& in these ranges affects 
the measurement by ±2.7% (±1.8%) for px > 4 (6) GeV. Varying the other parameters affecting 
the b-hadron and c-hadron sample compositions by their uncertainties has a smaller effect for 
both p T > 4 GeV and p T > 6 GeV (±0.7%, ±0.8%, respectively). 

b-quark properties: uncertainties in the production of B hadrons from the fragmentation of a b 
quark affect both the shape of the d xy distribution and the efficiency estimate. The systematic 
uncertainty is computed as the difference between the default result and those obtained with 
two different hadronization models in the PYTHIA simulation: the Lund symmetric |fTZ| and the 
Peterson [33 1 functions. Taking into account the effects on the b templates and those connected 
with the extraction of the efficiency, overall uncertainties of ±3.3% (px > 4 GeV) and ±3.6% 
(px > 6 GeV) are obtained. Using different PDFs to describe b-quark production in pp collisions 
has an effect of of ±0.9% (p T > 4 GeV) and ±0.5% (p T > 6 GeV). 

Light-meson decays in flight: muons from n and K decays have different d xy distributions. The 
shape is also different for muons from light mesons produced in the hadronization of a light 
quark, or from the decay of a heavy hadron. Given the uncertainties on the pion and kaon 
fractions in the simulation, we vary the relative amounts by ±30% and find a negligible effect 
on the final results. Similarly, we change the ratio of light mesons from heavy-flavor and light- 
flavor decays by ±50%, and observe a 2.5% (2.6%) change in the results for px > 4 (6) GeV. The 
generator-level filter applied to the simulated sample, requiring two muons to be produced 
within the tracker volume in each event, affects the shape and composition of the decays-in- 
flight template. The impact of the filter on the BB fraction is estimated by extracting the decays- 
in-flight template from an unbiased simulated sample in which only one generated muon is 
required to pass the filter and the other muon is used in determining the template. Repeating 
the analysis with this new template results in a 0.5% variation of the final result for both px 
selections. 

The total model-dependent systematic uncertainty, found by adding in quadrature the contri- 
butions listed above, is 5.5% for px > 4 GeV and 5.1% for px > 6 GeV. 
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7.2 Uncertainties on the impact parameter resolution 

The systematic uncertainty from the impact parameter resolution is determined by comparing 
the d xy distribution from prompt Y(1S) — > }i + u~ decay candidates reconstructed in collision 
data to the predicted distribution from MC simulation. A slight (p dependence in the determi- 
nation of the signed impact parameter with respect to the beam spot, where (p is the azimuthal 
angle of the muon track, due to the CMS tracker not being perfectly centered around the beam 
pipe, is not reproduced by the simulation. The combined effect from the misalignment and 
the different resolution in data and simulation is evaluated by an additional smearing of the 
impact parameter consistent with the observed differences between data and simulation. A 
further check to avoid the region dominated by the resolution has been performed by moving 
the lower bound of the fit range from to 40 pm. The maximum deviation from the default 
result found with these two methods is 2.7% (4.0%) for pj > 4 (6) GeV, which is taken as the 
systematic uncertainty due to the detector resolution. 

7.3 Uncertainties related to the Monte Carlo precision and the fit method 

There are four general categories of systematic uncertainty caused by the MC statistical preci- 
sion and the fitting procedure. These include: 

Monte Carlo precision: the likelihood fit is validated using a set of 500 parameterized simulated 
datasets, each with the same number of events as the data sample. The fit results from these 
datasets reproduce the input values with uncertainties consistent with those obtained in data, 
and the pull distribution is well described by a normal function. The r.m.s. of the results ob- 
tained for the BB fraction is 0.3% (0.7%) for pj > 4 (6) GeV, which are taken as the systematic 
uncertainties related to the finite simulated sample. 

Template parameterization: the d xy distributions in the simulated data used for the fit are smoothed 
using a superposition of a Gaussian plus one or two exponential functions, depending on 
the extent of the tail. The associated systematic uncertainty, evaluated by using different 
parametrizations, is equal to ±0.7% for both pj selections. The systematic uncertainty from 
the use of symmetrized templates is estimated to be ±0.6% (±0.7%) for pj > 4 (6) GeV, by 
comparing the results obtained in the simulation when a sum of symmetrized templates is 
used as pseudo-data instead of the usual randomized distribution. 

Bin size and fit upper bound: varying the d xy bin size in the range 0.002 — 0.008 cm accounts for a 
systematic uncertainty of 1.0% (2.1%) for pj > 4 (6) GeV, while varying the fit upper bound in 
the range 0.15 - 0.25 cm accounts for 0.3% (0.4%). 

Fit constraints: the BC, BD, and CD fractions are constrained in the fit so that their ratios with 
respect to the fitted BB fraction in the BC and BD cases, and the fitted CC fraction in the CD 
case, agree with the predicted values from the MC simulation, as described in Section [5J The 
uncertainties from this procedure include those on the rate of cc production from gluon split- 
ting and the muon misidentification rates in the simulation. To estimate the uncertainty, we 
vary the constraints on the fractions by ±50% around the simulation values, which induces a 
difference of 1.6% (1.2%) for p T > 4 (6) GeV in the fitted BB fraction. Since the 2D fit neglects 
the mixing of prompt and non-prompt muon components (PB, PC, PD), an additional system- 
atic uncertainty is computed by assigning to the BB fraction an uncertainty equal to the missing 
contributions, as found in the simulation, of 0.7% (0.6%) for pj > 4 (6) GeV. 

The total systematic uncertainty related to the fit method is found by adding the contributions 
in quadrature, which gives 2.2% (2.7%) for p T > 4 (6) GeV. 



7.4 Efficiencies from data and the dimuon invariant mass extrapolation 
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As a consistency check, an unconstrained ID fit is performed on the d xy distribution of the 



muons selected for the analysis, using the templates derived in Section 4.1 The results are in 
agreement within the quoted systematic uncertainty with those from the 2D fit. 

7.4 Efficiencies from data and the dimuon invariant mass extrapolation 

The statistical uncertainties of the efficiencies found from the T&P method of 6.0% (5.2%) for 
pi > 4 (6) GeV are taken as the systematic uncertainty on this procedure. 

The dimuon invariant mass distribution predicted from the MC simulation, scaled to the fitted 
fractions in the data, does not agree with the observed distribution within the uncertainties. 
Attributing the entire difference as being due to extra bb signal events, gives us the largest 
systematic uncertainty from this source of 1.1% (3.3%) for pj > 4 (6) GeV. 

7.5 Overall systematic uncertainty 

All the systematic uncertainties described so far are summarized in Table [5]and sum in quadra- 
ture to 8.9% (9.4%) for pj > 4 (6) GeV, with the larger contribution coming from the data-driven 
efficiency determination with the T&P method. The last source of systematic uncertainty to be 
considered is related to the integrated luminosity of the dimuon data sample, which is de- 
termined with a 4% uncertainty [16]. The total systematic uncertainty is therefore 9.8% for 
p T > 4 GeV and 10.2% for p T > 6 GeV. 

Table 5: Systematic uncertainties on the cross-section measurements in percent for the two pj 
limits. 



Source 


Uncertainty 




p T > 4 GeV 


pj > 6 GeV 


Model dependency 

Impact parameter resolution 

Monte Carlo precision and fit method 

Efficiencies and acceptance 


5.5 
2.7 
2.2 
6.1 


5.1 
4.0 
2.7 
6.2 


Total 


8.9 


9.4 



8 Results and comparison with QCD predictions 

The pp — > bbX — > jipX 1 cross section within the accepted kinematic range is determined from 
the observed number of dimuon events passing the event selection Nuy. t the fraction of signal 
events in the dimuon sample Jbb, the average efficiency for the trigger, muon identification, 
and event selection e, weighted by the pi and rj distributions, and the integrated luminosity C 
according to the relation: 



cr(pp -4bbX^ }i}iX',p T > 4 or 6 GeV, \tj \ < 2.1) 



■A/^ • /bb 
e-C ' 



By applying Eq. Q we measure: 

c-(pp -4 bbX -4 \i]iX!,pi > 4GeV, \rj\ < 2.1) = 

26.4 ± 0.1 (stat.) ± 2.4 (syst.) ± 1.1 (lumi.) nb 

and 

cr(pp -4 bbX -> wOC,pj > 6 GeV, \t]\ < 2.1) = 

5.12 ± 0.03 (stat.) ± 0.48 (syst.) ± 0.20 (lumi.) nb. 



(4) 



(5) 



(6) 
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The cross sections predicted by the leading-order pythia simulation are 48.2 nb for pj > 4 GeV 
and 9.2 nb for pj > 6 GeV, where the statistical uncertainties are negligible. That pythia pre- 
dicts a cross section value higher than the one measured in data has been seen in previous 
analyses [11], and is confirmed by our present findings. 

The next-to-leading-order event generator MC@NLO |34| is used to estimate the NLO QCD 
prediction for this measurement, with the CTEQ6.6 PDF and a b-quark mass of 4.75 GeV. The 
generator is interfaced with HERWIG [35 1 for parton showering, hadronization, and decays. The 
systematic uncertainty for this prediction is obtained by varying the b-quark mass between 
4.5 GeV and 5 GeV, and by changing the PDF to the MSTW2008 [36] set. The scale uncertainty 
is estimated by varying the QCD renormalization and factorization scales independently from 
half to twice their default values, as in Ref. fl37l . 

The predicted cross sections are: 

<7-mc@nlo(pP -> bbX -> fipX', p T > 4 GeV, \rj\ < 2.1) = 19.7 ± 0.3 (stat.) tH (syst.) nb (7) 
and 

^mc@nlo(pP -> bbX -> miX',p T > 6GeV, \ij\ < 2.1) = 4.40 ±0.14 (stat.) +J;J° (syst.) nb. (8) 

Both predictions are compatible with our results within the uncertainties of the NLO calcula- 
tions and the measurements. 

9 Summary 

A measurement of the inclusive cross section for the process pp — > bbX — > fijiX' at y/s = 7 TeV 
has been presented, based on an integrated luminosity of 27.9 ± 1.1 pb _1 collected by the CMS 
experiment at the LHC. Selecting pairs of muons each with pseudorapidity \rj\ < 2.1, the 
value cr(pp -> bbX -4 jtfiX') = 26.4 ± 0.1 (stat.) ± 2.4 (syst.) ± 1.1 (lumi.) nb was obtained 
for muons with transverse momentum pj > 4 GeV, and 5.12 ± 0.03 (stat.) ± 0.48 (syst.) ± 
0.20 (lumi.) nb for muons with pj > 6 GeV. This result is the most precise measurement of this 
quantity yet made at the LHC. 
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